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ABSTRACT 

We present exact solutions to the conformal Weyl gravity cosmological equations valid for 
both the matter and radiation dominated eras. We show that the cosmological models derived 
from this theory are not likely to reproduce the observational properties of our Universe. They 
fail to fulfill simultaneously the observational constraints on present cosmological parameters 
and on primordial light element abundances. 

Subject headings: cosmology: theory - early universe - gravitation - nucleosynthesis 

1. Introduction 

Conformal Weyl gravitational theory is based on the assumption of local conformal invariance of the 
space-time geometry. It has been proposed as an alternative gravitation theory (c.f. Adler 1982; Zee 1983), 
for the very high energy limit. Many attemps have been made at trying to break the conformal invarince 
dynamically so as to obtain Einstein gravity as a phenomenological theory. Unfortunately, no attempt has 
been successful as yet. Quantum fluctuations bring back the Hilbert-Einstein action with a repulsive rather 
than atractive gravity, when breaking the conformal symmetry dynamically. 

An alternative aproach to the symmetry breaking process has been recently proposed. Kazanas & 
Mannheim, in a series of papers, have suggested that Weyl gravity should be considered a theory in its 
own right, instead of looking for a low energy Einstein limit. They have found exact vacuum solutions 
(Mannheim & Kazanas 1989; 1990) and more recently, exact cosmological models, valid only for a radiation 
dominated universe (Mannheim 1992). 

Two important properties of this theory were found: both the cosmological constant problem and the 
flatness problem are solved naturally in the context of conformal gravity (Mannheim 1990; 1992) and the 
Newtonian limit of the theory gives rise to linear terms that could explain the flat rotation curves of galaxies 
without the assumption of a massive dark matter halo (Mannheim & Kazanas 1989). They also showed 
that the Schwarzschild solution can be obtained at length scales where the linear terms are unimportant 
(i.e. <C galaxy size), therefore fulfilling the same weak field observational tests that Einstein's theory does. 

Conformal Weyl gravity can be then considered as a valid theory of gravitation, at least in the weak 
field approximation. It is nevertheless necessary for the theory to be realistic to study the cosmological 
models and confront their predictions with the standard big bang cosmology. Mannheim (1992) has derived 
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cosmological field equations for the Robertson- Walker (RW) metric. He was able to obtain analytical 
solutions for a radiation dominated universe. The main result is that all possible models have an hyperbolic 
geometry (k < 0) but, contrary to the Friedman-Robertson- Walker (FRW) models, with oscillatory and 
singularity-free dynamics. Exact analytical solutions describing the whole thermal history of the Universe 
were not quoted. Besides, early universe nucleosynthesis in these models has not been studied as yet, 
although it is well known that Primordial Nucleosynthesis (PN) is one of the best ways to test the 
viability of alternative metric theories of gravitation (Will 1981; Dommguez-Tenreiro & Yepes 1987; Serna, 
Dominguez-Tenreiro & Yepes 1992). 

In this letter we present an exact analytical solution of the Conformal Weyl gravity cosmological field 
equations. Mannheim's solution for the radiation dominated era is just an approximation to this general 
one. We study those models that are compatible with present cosmological parameter values and compute 
the PN of these models. A comparison with the estimated primordial abundances of light elements shows 
that the models do not pass the PN test. 

In what follows we shall use units such that c = Ti = ks = 1. 



2. Cosmological equations 

Conformal Weyl gravity is described by the action 

I W =-a [ Cx^nC^y/^^X (1) 



where Ca^k is the Weyl tensor and a a dimensionless coefficient. The action for the matter content of the 
theory, Im, can be of the form (Mannheim 1992) 

Im = - j [S' I S,/2 + AS 4 - S 2 R»/12 + i?7"(*)(^ + M*))^ - hS^]^d A x, (2) 

where tp(x) is a fermion field representing all matter and S(x) is a conformal scalar field that yields 
symmetry breaking and renders the particles massive. R£ is the scalar curvature, T^x) the fermion spin 
connection and h and A are dimensionless coupling constants. 

Due to the conformal invariance of the theory the energy momentum tensor is kinematically 
traceless. Assuming a Robertson- Walker metric and the tracelessness of the following set of field 
equations is derived (Mannheim 1992) 

p+ \S 4 + S 2 (R 2 + k)/2R 2 = (3) 
3p- p-A\S 4 - S 2 (RR+ R 2 + k)/R 2 = (4) 

with p and p being the total density and pressure of a perfect fluid and R the expansion factor of the RW 
metric. Conformal invariance implies that S(x) = S is a spacetime independent constant. In addition, we 
also have the energy conservation equation for a perfect fluid 

3P + 3P+R dk~ = (5) 



3. Exact solutions 
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From the above set of equations we have been able to obtain exact solutions for the dynamical variables 
as a function of the temperature. From equations (3) and (4) and making use of the Hubble parameter 
H = R/R we get 

H =Y2 +Jr-2AS 2 -tf 2 (6) 

Introducing the definition of the decceleration parameter, q = —RR/R 2 = — (H/H 2 + 1), in the above 
equation 

2AS 2 = + qH 2 (7) 

This is a spacetime independent equation, therefore we can obtain the value of \S 2 evaluating eq. (7) at 
the current epoch. 

2AS 2 = ^+^ + g 5 2 (8) 



In what follows, a "0" subscript will stand for current parameter values. We are thus left with only one free 
parameter, the value of S. 

During most of the thermal evolution of the Universe (except in the particle-antiparticle annihilation 
era) the comoving entropy is conserved (RT =cte). In this case, eq (5) leads to (e.g. Weinberg 1972) 

f = ~ HT < 9 » 

that, in turn, can be used to transform the dependence of eq. (6) from time to temperature. Furthermore 
we can put p = p M + p R , where p M = pff (T/To) 3 is the matter density and p R = p R (T/To) 4 is the density 
of the relativistic fluid component. Also, p = p M + p R , with p M = and p R = p R /3. 

Changing to dimensionless variables % = H / Ho and T = T/To, and making use of eq. (7) and (9), we 
obtain 



d(lnft) g , 3fi J v 



M 



1 



1-T 3 + 2%(1-T 4 ) 
Po 



(10) 



d(lnT) H 2 8ttGS 2 H 2 

where we have introduced the density parameter i2ff = 8irGpff /3i? 2 . The solution of the above differential 
equation gives the temperature dependence of H. After a straightforward calculation we get, 

H = [-2BT 4 - 2AT 3 + (1 + q + 3A + 4B)T 2 -(q +A + 25)] 1/2 , (11) 

with A = 3i2ff /8irGS 2 and B = 3i2 R /8irGS 2 . It defines a set of cosmological models with only one free 
dimensionless parameter, GS' 2 . It must be noted that, contrary to the FRW cosmology, in this case qo is 
not directly related with S1q- 



From eqs. (3) and (4) we get: 



■ k Hp + p) 

X-R 2 "^ 2 —- ' (12) 



substituting H in the above equation, by means of eq. (6) and eq. (7) we obtain: 



k _ n 2 2p 3 Po + g0 2 



If we make use of the solution (11), after some arrangements we derive: 



R H T(l + g + 3A + 4B)i/2 ( 14 ) 



- 4 - 




Logio(T/T ) 

Fig. 1. — The Hubble parameter as a function of temperature for qo = 0.5, Slo = 1 and three different GS 2 
values: 10 5 (solid line), 10 10 (dotted-dashed line) and 10 18 (dotted line). For comparison, we also plot Ti(T) 
for the standard euclidean FRW model. 

which shows that only open (k < 0) universe models are possible since we know that Ho > 0. This result 
is directly related with the fact that the dimensionless constant A must be positive in order to have a well 
defined matter action (see. eq. (2)). From eq. (3) it can be seen that A > implies k < 0, since p is a 
positive quantity. 

There is not a general analytical solution for R(t). Nevetherless it can be easily seen that for a radiation 
dominated universe (A ~ 0) we can derive, from eq. (9) and (11), the same oscillatory solution found by 
Mannheim (1992). 

In our case, if we want to reach high enough temperatures for the PN to take place, we need that 
U 2 > 0, at least for T < 10 12 . This can be achieved only if A and B in eq. (11) are very small. This is 
equivalent to say that GS 2 must be very large. To show this behavior, we have plotted in Fig. 1 the exact 
solution Ti(T), given by eq. (11), for a wide range of GS 2 values and for the special case with qo = 0.5 and 
Slo = 1. The dependence of 7i(T) on the value of the cosmological parameters is negligible (if we restrict 
ourself to their observational limits) because the dynamical behavior is governed by the value of the scalar 
field S. 

From Fig. 1 we see that there always exists a maximum and minimum temperature corresponding to 
the minimum and maximum radius of the universe model given by eq. (14). 
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We can get the temperature at which H reaches the maximum value by finding the maxima of eq. 
(11). It is 



9A 2 1 
1 + ^ + — (l + go + 3A) 



-,1/2 



3A 
85 



(15) 



645 2 

For high GS 2 values, T Hmax ~ G5' 2 (S1^) _1 . In order to have temperatures > 10 9 K for the PN to take place, 
we should need GS 2 > 10 14 . For these values, % will be approximated by the asymptotic form 

Tiasym = + ?o)^ 2 ~ ?0 (16) 

We see the manifest linearity ofTi(T), contrary to the FRW model where % oc T 2 (see Fig 1). Thus, the 
ratio between the expansion factor of the standard FRW model and the one obtained here will be 

~Hfrw 



n 



oc T 



(17) 



This will be crucial for the viability of the models, since the PN abundances are very sensitive to large 
desviations from the Hfrw values at PN temperatures (Dominguez-Tenreiro & Yepes 1987). 

We have nevertheless computed the light element production in these models. We have used an 
updated version of the original Wagoner (1969) PN code that has been developped at Fermilab (Kawano 
1992). The code has been modified in order to account for the particular dynamics of the models described 
by this theory. To integrate the set of nuclear kinetic equations it is necessary to know how fast the universe 
cools. This is given by %(T). But the analytical solution (11) is not valid during the e + e~ annihilation 
phase, because the assumption RT = cte is no longer valid. It was necessary to numerically compute H(T) 
by eq. (13), using eq. (11) to set up the initial conditions at lower temperatures, where the analytical 
solution is valid. 



Table 1: Light element production in conformal Weyl gravity cosmological models. 



q fi /i 



n$h 2 



D/H 



(D+ 3 He)/H 



Y4 



He 



7 Li/ H 



0.5 



0.5 



0.25 



0.25 
2.5 
0.25 
2.5 



0.0075 
0.025 
0.25 
0.0075 
0.25 

0.5 

1.25 

2.5 
0.0075 

2.5 
0.0075 

2.5 



1.19 x 10 
3.13 x 10 
1.13 x 10 
1.42 x 10 
1.40 x 10 
x 10 
x 10 
x 10 
x 10 
x 10 
x 10 
5.24 x 10 



2.36 
5.69 
1.27 
7.08 
7.49 
6.38 



-20 
-20 
-19 
-20 
-19 
-19 
-19 
-18 
-21 
-19 
-21 
-19 



3.87 x 10" 8 

1.51 x 10" 8 
4.44 x 10" 18 
5.97 x 10" 8 
9.48 x 10" 18 
1.01 x 10" 17 
2.11 x 10" 17 

6.52 x 10" 17 
3.63 x 10" 8 
2.80 x 10" 17 
6.18 x 10" 8 
1.36 x 10" 17 



2.88 x 10" 
1.73 x 10" 
6.17 x 10" 
8.60 x 10" 
1.04 x 10" 
2.23 x 10" 
1.55 x 10" 
2.07 x 10" 
3.04 x 10" 3 
5.98 x 10" 2 

1.89 x 10" 3 
2.01 x 10" 1 



6.21 
1.53 

2.05 > 
2.85 
1.29 > 

4.06 > 

1.07 > 
2.01 > 
7.32 
4.74 > 
3.00 
5.31 > 



: lO" 9 

: lO" 8 
lO" 13 

: lO" 9 
10" 11 
lO" 13 
lO" 13 
lO" 12 

: lO" 9 
lO" 13 

: lO" 9 
lO" 14 



Note. — Abundances are by number relative to Hydrogen, except for 4 He that is given as mass fraction, h = Hq/100 km 
s" 1 Mpc" 1 ) 

In Table 1 we present the light element abundances for a wide range of baryonic content: from the 
lower observational bound Vt^h 2 = 0.0075 (Weinberg 1972) up to Vt^h 2 = 2.5 . This is equivalent to a 
baryon-to-photon ratio, r)s, ranging from r/s = 2.05 x 10 -10 up to 6.85 x 10 -8 . 
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There was no any combination of the parameter values that could give light element abundances within 
the observational limits (Boesgaard & Steigman 1985; Smith, Kawano & Malaney 1993). The PN process 
depends basically on the expansion rate and on the baryonic content. As we said before, the dynamical 
evolution of these models is not very much affected by the value of the cosmological parameters. Therefore, 
as can be seen in Table 1, the PN process is not very sensitive to the particular values of Ho, qo, or Slo • 

A correct 4 He abundance can be obtained in models with big amounts of baryonic matter 
(i2o h 2 ~ 0.5 — 2.5). The 7 Li abundance is also in agreement for a wide range of parameter values. But 
the main drawback of these models is the deuterium abundance. Because the universe described by these 
models expands much slower than in the standard FRW models, (see. eq. (17)), the nuclear reactions that 
burn the deuterium nuclei have more time to act consuming practically all of it. At the end of the PN 
process almost nothing of deuterium is left in contradiction with the current observational bounds that give 
a lower limit for the primordial abundance of D / H > 1.8 x 10 -5 (Smith, Kawano, & Malaney 1993). 

4. Conclusions 

We have obtained exact analytical solutions to the conformal Weyl gravity cosmological equations. 
We show that previous solutions, valid only for the radiation dominated era are only approximations to 
this solution. The main characteristic of the models are that they have an hyperbolic geometry but with 
an oscillatory behavior. Temperature is bounded between a maximum and a minimum value, with no 
singularity state. In order to have high enough temperature to allow PN to take place, we need to assume 
a very high value for the dimensionless parameter GS 2 . 

PN abundances in these models are incompatible with observational constraints for all range of 
parameter values. We conclude that conformal gravity cosmological models are very unlikely to give a 
realistic description of our Universe. 

We acknowledge M. Smith for kindly sending us a copy of the NUC123 PN code. We would also like 
to thank CP. Martin for very usefull discussions. 
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